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We apply the dynamical mean field theory (DMFT) in the iterative perturbation theory (IPT) to 
doubly degenerate e g bands and triply degenerate t2 S bands on a simple cubic lattice and calculate 
the spectrum and optical conductivity in arbitrary electron occupation. The spectrum simultane- 
ously shows the effects of multiplet structure and DMFT together with the electron ionization and 
affinity levels of different electron occupations, coherent peaks at the Fermi energy in the metallic 
phase and a gap at an integer filling of electrons for sufficiently large Coulomb U. We also calculate 
the critical value of the Coulomb U for degenerate orbitals. 
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I. INTRODUCTION 

Many physical properties in strongly correlated sys- 
tems have been studied extensively both experimentally 
and theoretically. These are not well described, how- 
ever, by the local density approximation (LDA) based on 
the density functional theory (DFT). LDA+U method 
has been developed beyond LDA, [1] and it works rea- 
sonably well for the insulating case near metal-insulator 
transition in strongly correlated systems. However, low- 
energy excitations in metallic phase near metal-insulator 
transition can not be treated by LDA+U method and 
more sophisticated approaches are needed. 

Recently the dynamical mean field theory (DMFT) 
has been developed and applied to some model sys- 
tems, which led us to a unified picture of both low- and 
high-energy excitations in the anomalous metallic phase 
near metal-insulator transition. [2] DMFT is based on 
the mapping of many body systems in bulk onto a single 
impurity problem, where the off-site Coulomb interac- 
tion can be neglected, subject to a self-consistency con- 
dition. The combination of DMFT with LDA is applica- 
ble to discuss realistic systems. [3, 4] The combination of 
DMFT with the GW approximation has been developed 
to include the off-site Coulomb interaction. [5] 

To solve the mapped single impurity problem within 
DMFT, one can use several computational schemes such 
as the quantum Monte Carlo method (QMC), [6, 7] the it- 
erative perturbation theory (IPT), [8, 9] the non-crossing 
approximation (NCA) [10] and the exact diagonaliza- 
tion (ED). [11] While QMC is not applicable in the low 
temperature limit, ED yields only a discrete number of 
poles for the density of states. Moreover, the exact cal- 
culations such as QMC and ED can only be suitable for 
simple model Hamiltonians. To carry out realistic calcu- 
lations in realistic materials, the approximate ones such 
as IPT or NCA would be more suitable due to its CPU- 
time, though IPT is not applicable for cases of large 
Coulomb interactions and NCA can not yield the Fermi 
liquid behavior at low energies and in low temperature 
limit. [12] 

The aim of the present paper is to apply the general- 



ized DMFT-IPT method for multi-bands on lattices [9] 
to the doubly degenerate e g bands and the triply de- 
generate t 2g bands on a simple cubic lattice and dis- 
cuss the fc-dependent spectra, the local Green's function 
and the optical conductivity near metal-insulator tran- 
sition. In Sec. II we will give the general formulation 
of the dynamical mean field theory and the optical con- 
ductivity. In Sec. Ill we will show our numerical results 
such as fc-dependent spectrum, the imaginary part of the 
local Green's function and the optical conductivity for 
both n = 0.5 and n < 0.5 cases. Finally we summarize 
our results in Sec. IV. 



II. FORMULATION 
A. Hamiltonian 

We proceed with the Hubbard-type Hamiltonian: 

e d c ]nuj c jma + ^mm' cr C jma C j' 'm' <J 

jma jj'amm' 

jmm' j(7m>m' 

(1) 

The first and second terms are the single-particle tight- 
binding Hamiltonian. Note that index j in the sum of 
Eq. (1) is running for correlated sites, {m} for orbital 
indices and {a} for spins. U is the on-site Coulomb in- 
teraction and the exchange interaction is neglected here. 

B. Dynamical mean field theory 

Lattice Green's function and local Green's function are 
defined as 

Grnam' a' (^i i^n) 

= [(iu n + M - eS)1 - M*0 ~ n^ n )] m \ m , a ,, (2) 
l r 

Gmtrm' a' {i^n) = 77 / G mam ' a ' (k, iiij n ) , (3) 
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where h(k), /i and 1 are the single-particle part of the 
Hamiltonian in fc-space, the chemical potential and the 
unit matrix. Here we neglect the fc-dependence of the 
self-energy E within the framework of DMFT. The chem- 
ical potential fi is determined to satisfy the Luttinger's 
theorem. [14] 

The IPT method was developed by Kajueter and 
Kotliar for non-degenerate orbital [8] and generalized by 
Fujiwara et al. for multi-orbital bands [9] on arbitrary 
electron occupation. In the IPT scheme, we assume the 
self-energy as 



X(iui n )=Un(N des - 1) + - y 

1 - B(iu} n )W>(iw n ) 



(4) 



where T.^(iLu n ) is the second-order self-energy, rid — 



= nNdcg is the total d-electron occupation 



number, n is the occupation number per each orbital and 
iVdeg is the degree of the degeneracy with respect to spins 
and orbitals. The coefficients A and B(iuj n ) are deter- 
mined by requiring the self-energy Yi(iu) n ) to be exact in 
the high-frequency limit (iuj n — ► oo), and in the atomic 
limit (U — ► oo) [9]. Y,^(iuj n ) is explicitly proportional 
to (iV dog - 1) and U 2 . 

The local Green's function is calculated in a form 

= / dkU ma (k)G a (k,iw n )U-^,(k), (5) 

a 
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G a {k,iu n ) = - 



iui n + n - e d - e a (k) - Y,(iu n ) ' 
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where the Hamiltonian matrix h(k) is diagonalized 
to be e a (k) by the unitary matrix U(k). The k- 
integration is carried out by using a generalized tetra- 
hedron method. [9] The 8 x 243 fc-points are used in 
the fe-integration within the whole Brillouin zone. For 
a paramagnetic and orbital degenerate system, the lo- 
cal Green's function becomes diagonal after the fc- 
integration. The Pade approximation is adopted for 
analytic continuation of the Green's function from the 
Matsubara frequencies iuu n to the real w-axis. 



From Eqs.(5) and (7), the optical conductivity is then 
calculated as 



Re 



de 



f{e)-f{e + u) 



ko act' 

xTr K(fc).^(fc)A a ,(fc). V (fc)], (9) 

where A a (k) = tf _1 (*0> I aa '{k,e u e 2 ) = 

D a (k,e 1 )D a >(k,e 2 ), E™ m ' = 8 mm ,5 ma , andD Q (fc,w) = 
— -ImG a (fc, co). The fc-integration in Eq. (9) should be 
carefully carried out. 

The fc-integration is able to be carried out simply 
by the generalized tetrahedron method if the integrand 
has one pole. Since the two-body Green's function 
I aa '(k, ei, e 2 ) has two poles, we should rewrite it to sep- 
arate into two partial fractions as 

I aa > (fe, ei, e 2 )=^Re [G a (k, ei)G*,(fc, e 2 ) 
-G a (k,e 1 )G a/ (k,e 2 )} 

= i Re ( G Q (fc, ei )-G;,(fc, £2 ) 

2 Ua 1 (fc,e 1 )-{G;,(fc,e 2 )}-i 



G a (fc,6i) - G a /(fc,e 2 ) 
Ga 1 {k 1 e l )-G-J{k 1 e 2 ) 



(10) 



In case a = a' , the denominator of Eq. (10) is actually 
fe-independent and the generalized tetrahedron method 
can be applied. 

When lj <~ 0, we should treat Eq. (10) separately for 
a = a' and a ^ a' . When u ~ and a = a' , Eq. (10) 
is then rewritten as 



I aa (k,e, e)=-Re 



Z(e)= 



\G a (k,e)\ 2 + Z(e) 



dG a (k,e) 



de 







duj 





(11) 

(12) 



When oj ~ and a 7^ a', Eq. (10) should be rewritten as 



C. Optical conductivity 

From the linear response theory, the optical conduc- 
tivity is given as 



Re o- MM / (w) = J dt 



/(e) -/(e + w) 



xTr [Im G(k, e)j M (fc)Im G(k, e + (fe)] (7) 
The matrix element of the current operator is defined as 

aimraVu 

3™ (k)=-e dk 1 > {n = x,y,z). (8) 



fact' (^j ^7 ^ 

1 



Re [{(e - ReE(e) - e a (k))(e - ReE(e) - e a ,{k)) 



+ (Im E(e)) 2 + i(e a (fe) - e a > (fc))Im E(e)}- 1 
G a (fc,e)G a ,(fc,e) 
e a (fe) - e a '(fe) 



(13) 



In the above two cases, the resultant I aa >(k, e, e) can be 
calculated by using the generalized tetrahedron method. 

For cases of t 2g bands, the fe-dependent Green's func- 
tion is diagonal and this is also true for I aa > (fe, e, e) . How- 
ever, for for e g bands, the fc-dependent Green's function 
has non-zero off-diagonal elements and also I aa /(k, e, e) 
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FIG. 1: The imaginary part of the fc-dependent Green's function G(k, uj) (shaded regions) of e g bands and t2 S bands on a simple 
cubic lattice. The energy bands e(fc) shifted by ±17/2 at high-symmetry lines are also shown. The high-symmetry fc-points 
are M(f , f ,0), X(f ,0,0), r(0,0,0), R(f ,f ,f ). (a) e g bands, U = 3.0. (b) e g bands, U = 2.0. (c) e g bands, [7 = 1.0. (d) t 2g 
bands, U — 3.0. (e) fa s bands, U — 2.0. (f) t2 S bands, U = 1.0. The inverse temperature /3 — 30. In the e g bands 
cases (a)~(c), red and green shades represent Gn(fc,w) and G22(k,u>) of the diagonalized fc-dependent Green's function 
matrix, respectively. In the t2 g bands cases (d)~(f), red, green and blue shades represent G VZj y Z (k,u>), G Z x,zx(k,ui) and 
G xy ,xy(k,Lo), respectively. 



has. Then the denominators of Eq. (10) in case of 
a 7^ a! for e g bands are fc-dependent and, therefore, 
fc-integration is not able to be carried out by using 
the generalized tetrahedron method. For this case, fc- 
integration should be carried out as the sum of the inte- 
grand at each u> + i5. 



D. The energy bands 



We have applied the above described calculation 
scheme to the doubly degenerate e g bands (iVd eg = 4) 
and the triply degenerate t2 g bands (iVdeg = 6). We con- 
sider Slatcr-Koster-type tight-binding Hamiltonian on a 
simple cubic lattice. The effective hopping integrals are 
assumed only between the nearest neighbor pairs and 
V dd a = 1/3, ViAr = -2/3V dda and V ddS = l/6V dda . The 
relationship among Vdda-, V ddTr and V dd s is the scaling 
properties of bare two-center integrals in the LMTO 
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FIG. 2: The imaginary part of the local Green's functions in a half-filled case n = 0.5. (a) Doubly degenerate e g bands of the 
half band width D — 7/6. (b) Triply degenerate t2 S bands of the half band width D = 1. The inverse temperature f3 = 30. 



method. Here the half band width of e g orbitals and t2 g 
orbitals are D &s = 3{V dda + V ddS ) = 7/6 and D t2s = 
2(V dd7T — V dd s) = 1, respectively. The off-diagonal el- 
ements of h{k) of arbitrary k are non zero for the e g 
bands. On the other hand, those for the t 2g bands are 
always zero. 

III. NUMERICAL RESULTS 
A. fc-dependent spectrum 

The imaginary part of fc-dependent Green's function 
G(k,u>) is shown in Fig. 1 (shaded regions) with the en- 
ergy bands e(fc) (full lines), for f3 = 30 and the electron 
occupation n = 1/2 (a half-filled case). The energy ze- 
roth is set to be at the chemical potential /i. Energy 
bands e(fc) are shifted by ±{7/2, corresponding to lower 
and upper Hubbard bands. 

Within the Fermi liquid theory, the fc-dependent 
Green's function behaves as 

G(fc, u) ~ Z/{ul - Zh(k)} (14) 

near uj ~ in a metallic case, where Z = (1 — 1 is 
a renormalization factor. Then the width of the coher- 
ent peak at the Fermi energy is scaled by 1/Z. These 
values become constant within the DMFT, since the fc- 
dependence of the self-energy E is fc-independent. 

In the e g band case, peaks appear at the Fermi en- 
ergy for U = 1.0 and U = 2.0, which shows that the 
system is in a metallic phase. For U — 1.0, those peaks 
strongly depend on fc and the spectrum becomes sharper 
rapidly near the Fermi energy. On the contrary, for 
U = 2.0, those peaks hardly depend on fc. With in- 
creasing Coulomb interaction in a metallic phase (U = 



2.0), coherent peaks appear. At U — 3.0, the spectrum 
splits into upper and lower Hubbard bands. This results 
from the combination of electronic structure calculation 
and DMFT. The model calculation of DMFT with only 
assuming the density of states is not able to give the fc- 
dependent spectra. 

It must be noted that, in t2 g case, the spectrum shows 
already the coherent peak hardly depending upon fc at 
U = 1.0 and it splits into upper and lower Hubbard bands 
at U = 2.0. 



B. Imaginary part of the local Green's function 

1. n = 1/2 case 

The imaginary part of the local Green's function 
G mm(T (w) is shown in Fig. 2, for (3 = 30, with the electron 
occupation n = 1/2 (a half- filled case). The positive and 
negative energy regions, corresponding to the spectrum 
of excitation from occupied to unoccupied states (elec- 
tron ionization spectrum) and that from unoccupied to 
occupied states (electron affinity spectrum). 

In a half-filled case, the spectrum consists of the upper 
and lower Hubbard bands with the electron-hole symme- 
try. At the critical region (here, U ~ 2.0 for e g bands 
and U ~ 1.0 for ti g bands) one observes a sharp coherent 
peak at lo ~ 0. In an insulating case, the self energy 
shows the singular behavior E ~ — near ui ~ 0. On the 
contrary, In a metallic case, the self energy behaves as 
£ — Uji(N deg — 1) ~ uj near lo ~ 0, which is typical for 
the Fermi liquid. 

As seen in Fig. 2, the critical U/D ratio is different 
between e g bands and t 2g bands; (U/D) cr ~ 1.71 for 
e g bands and (U/D) cr ~ 1.0 for t 2g bands. Figure 2 
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FIG. 3: The electron occupation dependence of the imaginary part of local Green's functions for the triply degenerate £2 g bands 
of the half band width D — 1 and the inverse temperature /3 = 30. The numbers in the figure denote the electron occupation 
per each orbital. Red line shows the case of the (almost) integer filling of the total d-electron number n = 1/6 (rid = 1), n = 
1/3 (n d = 2) and n = 1/2 (n d = 3). 



also shows that in the insulating case, the band width 
is shrinking more at the same U/D ratio for t2 g bands 
than that is for e g bands. This is because that the self- 
energy is much larger for t2 S bands than for e g bands 
for the same U/D value due to the degree of the de- 
generacy with respect to spins and orbitals, N^eg- The 
second-order self-energy is proportional to (Nd eg — l)U 2 
and it contributes 5/3 times larger for t2 g bands than 
for e g bands. This fact means that the critical U/D ra- 
tio of metal-insulator transition decreases with increasing 
degeneracy, which contradicts to the conventional under- 
standings that the critical U / D ratio would be the same 
as in the s-orbital case, [15] or that the critical U/D ra- 
tio would increase with it. [16] These contradiction could 
originate from the Hamiltonian on a lattice or the sym- 
metry type of the orbitals. 



2. n < 1/2 case 

Figure 3 shows the electron occupation dependence 
of the imaginary part of the local Green's functions for 
t2 g bands with (3 = 30. (See Ref. 9 for e g -band.) In 
case of n — 1/6 (rid = 1), n = 1/3 (n<j = 2) and 



n = 1/2 (rid = 3), one observes the coherent peak at 
uj ~ for U = 1.0, which is the critical region of metal- 
insulator transition. We also observes the gap for the 
same n's at ui ~ for U = 2.0 and system becomes in- 
sulator. These occupations correspond to the cases of 
the integer filling of the total d-electron number. At 
the integer filling case of rid — i (idnteger), the states of 
rid = i is completely filled and the states of rid = i + 1 
is completely empty. Then the gap opens between the 
ionization and affinity spectra. 

With small hole doping S, in which electron occupa- 
tion rid = i — 5, the spectrum changes drastically at 
U = 1.0 and U = 2.0. Since the ground state is the 
mixture of the states of rid = i and rid — i — 1, a sharp 
peak appears below the lower Hubbard band. This 
satellite peak corresponds to the electron ionization of 
rid = i — * i — 1. The ionization and affinity levels for 
both rid — i and rid = i — 1 states appear in the spec- 
trum. These satellite structures of all configurations of 
rid = 0, 1, 2, ... appear but their intensity depends on the 
occupation number. With increasing hole doping i5, the 
electron ionization spectrum gradually increases. On the 
contrary, the electron affinity spectrum loses the inten- 
sity. When the occupation approaches rid — i — 1 from 
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FIG. 4: The optical conductivity in a half-filled case n = 0.5. (a) Doubly degenerate e g bands of the half band width 
D = 7/6. (b) Triply degenerate <2 g bands of the half band width D = 1. The inverse temperature j3 = 30. The arrows 
indicate the height of the Drude peak. 



above, the electron affinity spectrum n d = i —> i + 1 dis- 
appears. 

From above discussion, we have had success of the 
spectrum calculation which simultaneously shows the ef- 
fects of multiplet structure and DMFT such as the elec- 
tron ionization and affinity levels of different electron oc- 
cupations, coherent peaks at the Fermi energy in the 
metallic phase and a gap at an integer filling of electrons 
for sufficiently large Coulomb U . 



C. Optical conductivity 

1. n = 1/2 case 

Figure 4 shows the optical conductivity in a half-filled 
case. In a metallic case, (U = 1.0, 2.0 for e g bands, U — 
0.5, 1.0 for £ 2g bands) three peaks appear: (1) the Drude 
peak at ui ~ 0, (2) a peak of the width 2D at ui < U and 
(3) a peak of the width AD at uo = U. The peak (2) 
comes from the transition from the coherent state to the 
upper Hubbard band or that from the lower Hubbard 
band to the coherent state. The peak (3) at U comes 
from the transition from the lower Hubbard band to the 
upper Hubbard band. The peaks (1) and (2) result from 
the DMFT which includes the coherent state. 

With increasing Coulomb U in metallic phase, con- 
tributions of the peaks (1) and (2) decreases. In the 
insulating case, (U = 3.0 at e g bands, U = 2.0, 3.0 at £ 2g 
bands) the peaks (1) and (2) disappear and only that of 
(3) remains. Figure 4 is consistent with the fact that the 
critical U/D ratio is different between e g bands and 
bands [(U/D) cr = 1.71 for e g bands and (U/D) cr = 1.0 



for t2 g bands]. 



2. n < 1/2 case 

Figures 5 and 6 show the optical conductivity in cases 
of arbitrary electron occupations. In case of the inte- 
ger filling of the total d-electron number, such as n = 
1/6 (n d = 1), n = 1/3 (n d = 2) and n = 1/2 (n d = 3), 
for U = 2.0, the Drude peak disappears and one observes 
only a peak at w = U, which comes from the transi- 
tion from the lower Hubbard band to the upper Hubbard 
band. 

With small hole doping S, in which electron occupation 
rid = i — 5 (i=an integer number), the Drude peak ap- 
pears and the peak at lo = U shrinks. One also observes a 
peak which is separated from Drude peak at u ~ 0.4 for 
both U — 1.0 and U = 2.0. With increasing hole dop- 
ing <5, this peak shifts toward lower energy and becomes 
sharper and merges into the Drude peak, because it 
comes from the transition from the lower Hubbard band 
to the coherent state. At this range, the hight of the 
Drude peak becomes higher and sharper. Also we can 
see that the intensity at the lower energy side of the peak 
at lo = U decreases simultaneously, and this is because 
that this side of the peak at u> = U were mainly con- 
tributed by the transition from the coherent state to the 
upper Hubbard band. When the occupation approaches 
rid = i — I from the above, the Drude peak disappears 
again for U — 2.0. 

Optical conductivity successfully shows the ionization 
and affinity transition and also metal-insulator transi- 
tion. This is, of course, consistent with the results of the 
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FIG. 5: The electron occupation dependence of the optical conductivity at U = 1.0 for the triply degenerate t2 g bands of the 
half band width D = 1. The inverse temperature /3 = 30. The numbers in the figure denote the electron occupation per each 
orbital, (a) < n < 1/6 (0 < n d < 1), (b) 1/6 < n < 1/3 (1 < n d < 2), (c) 1/3 < n < 1/2 (2 < n d < 3). The arrows indicate 
the height of the Drude peak. 




FIG. 6: The electron occupation dependence of the optical conductivity at U = 2.0 for the triply degenerate t2g bands of the 
half band width D = 1. The inverse temperature /3 = 30. The numbers in the figure denote the electron occupation per each 
orbital, (a) < n < 1/6 (0 < n d < 1), (b) 1/6 < n < 1/3 (1 < n d < 2), (c) 1/3 < n < 1/2 (2 < n d < 3). The arrows indicate 
the height of the Drude peak. 



spectrum. 

IV. CONCLUSION 

In conclusion, we have applied the generalized DMFT- 
IPT method to doubly degenerate e g bands and triply 



degenerate t2 g bands. The fc-dependent spectrum shows 
the resonance behavior near the Fermi energy at smaller 
U and the coherent peak near metal-insulator transi- 
tion. The spectrum shows that the critical U/D ra- 
tio of metal-insulator transition is different between e g 
bands and £2g bands and this originates from the degen- 
eracy Ndeg- In the arbitrary occupation case, the spec- 
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trum shows ionization and affinity levels, coherent peaks U. The optical conductivity gives a spectrum consistent 
at the Fermi energy in a metallic phase and a gap at an with the above discussions, 
integer filling of electrons for sufficiently large Coulomb 
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